On the nature of ferromagnetism in diluted magnetic 
semiconductors: GaAs:Mn, GaP:Mn. 

V.A. Ivanov 1,2 , P.M. Krstajic 1 , F.M. Peeters 1 , V. Fleurov 3 , and K. Kikoin 4 

1 Department Natuurkunde , Universiteit Antwerpen (VIA), 
Universitetsplein 1, B-2610 Antwerpen, Belgium 
2 N.S. Korsakow Institute of General and Inorganic 
Chemistry of the Russian Academy of Sciences, 
Leninski prospect 31, 117 907 Moscow, Russia 
3 School of Physics and Astronomy, Tel Aviv University, Ramat Aviv, 69 978 Tel Aviv, Israel 
4 Physics Department, Ben-Gurion University, 84 1 05 Beer-Sheva, Israel 

(Dated: February 1, 2008) 

A microscopic Hamiltonian for interacting manganese impurities in diluted mag- 
netic semiconductors (DMS) is derived. It is shown that in p-type III-V DMS the 
indirect exchange between Mn impurities has similarities with the Zener mechanism 
in transition metal oxides. Here the mobile holes and localized states near the top of 
the valence band play the role of unoccupied oxygen orbitals which induce ferromag- 
netism. The Curie temperature estimated from the proposed kinematic exchange 
agrees with recent experiments on GaAs:Mn. The model is also applicable to the 
GaP:Mn system. 

PACS numbers: 71.27.+a, 75.30.Hx, 75.50.Pp 

The discovery of ferromagnetism (FM) with Tc=110K in Gao.947Mno.o53As |IJ stimulated 
the systematic study of the III-V dilute magnetic semiconductors (DMS). Recently, above 
room temperature FM order was found in p- type GaP and GaN doped with Mn [g. 

Most of the existing theories of FM in III-V DMS are based on semi-phenomenological 
models, which postulate the existence of local magnetic moments on the Mn sites, indirect 
exchange between these moments and holes in the valence band of the host crystals (see, 
e.g. ||) and sometimes emphasizes the role of shallow acceptor levels [|J. In this paper 
we present a microscopic model, which takes into account the origin of localized magnetic 
moments and shallow acceptor levels induced by the Mn impurity and derive the effective 
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kinematic exchange from the generic two-impurity Hamiltonian. This mechanism cannot 
be reduced to any previously proposed models, but it is based on Zener's idea of double 
exchange via unoccupied p-orbitals. 

In the present paper the basis for a microscopic description of an isolated magnetic 
impurity is the Anderson Hamiltonian || which was modified in Refs. [0, [8] to the case 
of a semiconductor host. A two-impurity generalization of the Anderson model for metals 
was proposed in Ref. || (hereafter referred as AA). An indirect exchange between magnetic 
moments arises as a result of virtual electron transitions into unoccupied states shared by 
two impurities. 



It is known (see Refs. JL0], [TTJ]) that 3d impurities create both deep localized states in the 
energy gap and resonance states in the valence and conduction bands. Three types of states 



are generated by the 3d impurities in zinc blende semiconductors [12, 13, 14]: non-bonding 
states retaining the angular e-symmetry of the states in a cubic crystal field and two types of 
t2 states (bonding and antibonding). The latter arise due to a strong hybridization between 
the atomic t-i orbitals and the p-states of the same symmetry belonging mainly to the heavy 



hole (hh) band [T^, 11]. They are called crystal field resonance (CFR, predominantly d- 
type) and dangling bond hybrid (DBH, predominantly p-type). One of these states, as a 
rule, gives an impurity level in the energy gap and another manifests itself as a resonance 
within the valence band. According to its position in a series of transition metal elements, 
the Mn-impurity in the Ga-site should have the configuration 3d 4 . However, the Mn ion 
retains its fifth electron in the 3d shell because of a specific stability of a high-spin half-filled 
state d 5 and the impurity state is Mn(3d 5 p), p is the bound hole state. Actually, in GaAs:Mn 
the mobile hole concentration is about 30% of the nominal Mn concentration JT5[ . 

The p<i-hybridization together with the Anderson-Hubbard repulsion U is eventually the 
source of the magnetic interaction in DMS. To describe the indirect exchange, we start with 
the single impurity resonance scattering model 0, || for ti electrons and generalize it to the 
case of two impurities along the lines of the AA approach [[|. The 'passive' non-bonding 
e-states contribute to the localized moment, but do not participate in the indirect exchange. 
Therefore, the minimal two-impurity Hamiltonian involves t 2 -electrons: 

V 

where only hh states are retained in the band Hamiltonian. Here Cp ha (c p h a ) is the creation 



H = J2 4 c lh* c ph° + £ (V^cl^e^ + h.c) + £ (E d n° + \n°< ) , (1) 
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(annihilation) operator of a hh with the momentum p and spin a. The second term describes 
the resonant impurity scattering induced by the pd-hybridization V p d- The last term contains 
the atomic energy levels of the localized d electrons with the ^-electron occupation operator 
nl = d\ a di a of the Mn impurity in the Ga sites labelled by % — 1, 2. 

The system of Dyson equations for the d-electron Green's functions G dii , (i, i' = 1, 2) has 
the following form: 

G%g(e) = g ia (e) U, + V 2 ^ (e)G^(e)J , (2) 

where gi a {s) = {is — Ed — Un^ <7 )~ 1 is the single site Green's function for a t2 electron, 
Lij(s) = J2 P e ~*^ '( R '~ Rj )(2£: — £p) _1 is the lattice Green's function for hh. The momentum 
dependence of the hybridization matrix elements is neglected, i.e. V p d ~ V. The solution 
of the system of Eqs © is G a dii = [g£(e) - V 2 L^{e)] /R%e), G a dij = V 2 L^s)/R%s), 
(i = 1, 2; j = 2, 1) and the two-impurity levels are found from 

^(e) = II - V 2 L^)] ~ V'L^L^e) = 0. (3) 

t=l,2 

The zeros of the expression in the square brackets describe the impurity o?-level renormal- 
ized by their hybridization with the hh band. For sufficiently large hybridization (or narrow 
valence band) the DBH states arise above the top of the valence band, whereas the CFR 
levels appear deep in the valence band below the bottom of the hh subband (the left panel 
of Fig. |l|). The inter-site interaction V 4 \L1 2 \ 2 results in mixing of the CFR and the DBH 
states belonging to the two impurities (right panel of Fig. 0). 

The occupied CFR levels E ia correspond to the states d 5 /d 4 of the Mn ions, whereas 
the empty d 6 /d 5 CFR levels (E^) are shifted to the conduction band by the Anderson- 
Hubbard repulsion U (Fig. HD responsible for the spin-dependent inter-impurity interaction 
and eventually for the FM order. 

Both in GaAs:Mn and GaP:Mn the deep CFR states are completely occupied, and the 
DBH states in the energy gap are empty (see, e.g., Refs. |TT| , PIP- Two competing mech- 
anisms of the magnetic interaction arise, since the indirect exchange between neighboring 
impurities involves either the empty states near the top of the valence band and the empty 
(d e /d 5 ) CFR levels. To determine the type of magnetic ordering, we consider the indirect 
exchange between two neighboring magnetic ions and calculate its sign and magnitude. 



FIG. 1: Left panel: graphic solution of equation (|3|) for the bonding CFR and antibonding DBH 
levels. Right panel: energy levels in GaAs:Mn. The CFR d-levels d 5 /d 4 (denoted by R 1>2 ) of each 
impurity, lie below the hh band. The DBH levels (energies E\ a , e<io) are split from the hh band 
and form localized (acceptor) levels in the energy gap. The CFR levels d 6 /d 5 R-^-2 lie high in 
the conduction band. 

The impurity related correction to the energy of the system is given by the standard 
formula 

1 poo 1 

E ma gn = _ Im / e ATr G[e - z5sign(e - e F )]de - - ^ Un di ^n dil , (4) 

7T J-oo L T 

where G (E) = (E — H)^ 1 is the full Green's function. We estimate first the contribution of 
the CFR levels. In this case the two-site lattice Green's function L 12 (e) ~ L 12 (E ba ) = L is 
approximated by its value for the CFR d -level position 



E ia = E d + V 2 L u {E ia ). (5) 
The quantity L depends exponentially on the intersite distance R\2'. L ~ 
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exp(— KbRu)/ (KbRuw) where Kb = \j2m{e bh — E ia )/h, with Sbh standing for the bottom 
of the hh band, w is the hh bandwidth, E ia is the energy of the CFR level below the bottom 
of the hh band. 

Due to the on-site repulsion U, the structures of discrete CFR levels differ for FM and 
antiferromagnetic (AFM) orientations of neighboring impurity spins. For the AFM align- 
ment the occupation numbers are nif = n 2 i = 1, nil = n%<\ = . As a result, the tunneling 
processes which influence energy positions of the occupied states involve large U, and the 
secular equation @ yields the following expressions for the occupied impurity levels (see 
Fig. H): Ebi ~ Ei^ — J a, E b i ~ E 2 ± — J a, where the indirect exchange parameter 

J A = V 4 L 2 /U . (6) 

is the Anderson superexchange, which favors AFM order in transition metal oxides. 

In the FM case: = n 2 -\ = 1, nij = n 2 [ = 0. Then the states, which mediate exchange 
between the impurity spins, are the mobile hole states below the top of the valence band 
and the localized DBH states above its top. The role of these states in the DMS is the 
similar to the that of the empty atomic states in the conventional Zener double exchange 
mechanism for (La, A 2+ )Mn03 [[J. An important difference is that in Zener's case the Mn 
ions are in different valence states (Mn 3+ and Mn 4+ ). In other words, one of the two levels 
* s em Pty- Since in our case both of these states are occupied, the Zener mechanism in its 
original form does not work. Formally, one gets a pair of bonding/antibonding CFR states 
E(b, a )i = Erf ± Jz from flj|) without any energy gain, since both levels are occupied. It will 
be shown below that FM order arises only at a finite hole concentration in the valence band. 

Impurity related corrections to the band energy are obtained from the integrations in Eq. 
©■ ( c f- R e f- ID where a similar procedure was carried out for a single impurity). Then, 
the variation of the total band energy due to the two-impurity scattering is 

AE b = --Im r de\n{[g^{ie) - V 2 L lx {ie)f - V*Lj 2 (ts)}. (7) 

Here i-y(ze) = Pij(e) + |ry(e) and 

Pn(e) = P 2 2{e) = / duj , P 12 (e) = / duj- 

J e — oo J 



kR 12 e — u 

r 12 (e) w 2np(e) smkR 12 /kR 12 , r u (e) w np{e) = T 22 (s). (8) 

Only the spin-up (majority spin) band contributes to AEb- Here and below the spin in- 
dex is omitted for the sake of brevity. The value of the wave-vector k is found from the 
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equation e = £hh(k), where £hh{k) is the hh energy dispersion. Since V/ [g^ 1 — V 2 Pnj = 
V I {e — Ed — V 2 Pn (e)) <C 1 one can extract from Eq. ([?]) the spin-dependent hh band 
contribution to the exchange energy, which reads 



, ox f f° r 12 (e)P 12 (e) , „_ P l2 (e i )P{ 2 (e i 



where "0*1,2 are the vectors of Pauli matrices. The terms in curly brackets correspond to the 
contribution of mobile and localized holes. The factor 3 appeared due to the degeneracy of 
the localized acceptor p-levels £j. The function i^-(e) = dP i j(e)/de is negative at e = e$; (see 
left panel of Fig. [I]). 

Eq. (|^) is the main result of our theory from which we obtain Tc- In the evaluations of 
Anderson-type and Zener-type coupling constants we use the estimates: Py ~ w~ x , ~ 
e^w-V 2 , e f - E d - V 2 P u {e f ) = 4aV 2 /w with a < 1 (see left panel of Fig. [I]). Then 
one finds from Eqs. (j|) and (|) that J a ~ V A [exp (—HbRu) / (i^bRn)] 2 /(Uw 2 ), with «;& = 



2m(£ bh —E da )/h, J F ~ 2ef{ef/w) 1 I 2 / (Aa) 2 + Ep/ (4u>) and FM pairing is realized 
provided Jf > Ja- 

To calculate the exchange energy Eq.®, one needs the dependence of the Fermi level 
on the Mn concentration Ep(x). It is governed by the equation x s = 2 f® p p v (E)dE, where 
x s is the hole concentration per site. We assume a semi-elliptical density of states, p v (e) = 



2 [ttw ] y — e(e + w)8(—e), where w = 2.9 eV p!7| . The hole concentration per site, x s , is 
proportional to the hole density per volume ph- x s = l/8a 3 ph, as there are four III-V pairs 
in a unit cell volume a 3 in zinc-blende structures. Based on the experimental data Ph{x) of 
Ref.[F6 we used a polynomial fit for the hole density dependence on the Mn concentration 
in GaAs:Mn. 

To compare the values of J a and Jp we use the calculated Ef{x — 5.3%) = — 50meV, 
U ~ 4.5 eV, the hybridization parameter V = 1.27eV obtained from (§]) for the acceptor 
level Ei = 85 meV ( £. xp = 110 meV 0) and the CFR level E t = -3.0eV (E- xp = -3.4eV 
IH). At these values a = w [e F - E d - V 2 P 11 {e f )} / {AV 2 ) = -0.32 and the ratio r = 2.13 
justifies the dominance of FM coupling in GaAs:Mn: Jp > Ja 

Eq. (|) allows one to compute the Curie temperature T c (x) = E ex (x)/k B . The results 
are presented in Fig. [|. The calculated dependence Tciph) is given in Fig. ^| for different 
values of the Mn concentration. 



FIG. 2: The dependence of the kinematic exchange (left axis) and the Tq (right axis) on the 
impurity concentration. Experimental results (filled circles with error bars) are taken from Ref. 
[16]. The solid and dashed curve are obtained using the experimental obtained central and upper 



values of the error bars for the hole density, respectively of Ref. |T 



Although we neglected in this paper formation of an impurity band around x ~ 5% our 
results for Tq(x) and Tc{ph) are in good agreement both with available experimental data, 
and the theoretical phenomenological estimates. The dielectric structure parameters such 
as DBH, lower CFR energies and V are calculated self-consistently (Eq. fl5|) and see the 
left panel of Fig.[T[) for the graphical solution). Even our estimated value of V — 1.27 eV 
is in a reasonable agreement with the hybridization parameter ~ 1.4 eV which can be 
extracted from the data of Ref.[^[. This implies that the present theory, in fact, has no 
real fitting parameters. Note that in Ref. [f2l]] the authors took into account the magnetic 
correlation mediated by holes originating from shallow acceptors without a proper regard of 
the role of the impurity levels, which we consider to be important, and later calculated Tq 
by minimizing the Landau free energy. 



FIG. 3: The predicted dependence Tq {Ph)i at the fixed manganese concentration x = 5.3% and 



7%. (solid lines). The dashed line are the theoretical result of Ref. [20], derived from the 
phenomenological Landau free-energy functional within the framework of the Zener model. The 
dash-dot line is the recent result [24] obtained by using the exchange-enhanced mean-field theory. 



The point at x = 5.3% corresponds to the experimentally found [16] maximal value of Tq for 
GaAs:Mn. 

In conclusion, we proposed a microscopic model for a double exchange in p-type III- 
V:Mn DMS based on the known mechanism of interaction between substitutional transition 
metal impurities and the host semiconductor [Kj, [|TT|. A presence of holes is crucial for 



the FM double exchange between neighboring Mn ions. Our model there does not require 
an introduction a adhoc phenomenological pd- or RKKY-type exchanges. The source of the 
magnetic coupling is the energy gain in the kinetic energy of holes in case when neighboring 
ion spins are parallel. Among existing approaches the closest one to ours is the LSDA+U 
method|pl|, in which the system of t^cf- and ea- levels and the hole pocket for majority 
spin density of states correlates with our spectrum presented in Fig. [I]. However, instead of 
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extracting the j?d-exchange from this band structure we calculated the genuine Zener-like 
exchange which emerges since the band energy is lower in the FM case then in the AFM 

case 



22fl . In our model this exchange is defined by Eq. (|]). Similar ideas were applied to 
DMS CdGeP2:Mn (||23|), where the interplay between CFR and DBH states turned out to 
be crucial for the onset of FM order. The theory may be applied to GaP:Mn which has a 



similar structure of chemical bonds around the Mn ion |TJ| ; however, experimental data on 
the CFR and DBH states as well as the dependencies Ph (x) and Tq(x) are not available. 
The proposed kinematic exchange can also be applied to other DMS including GaN:Mn and 
A 2+ GeP2:Mn, which is left for future work. 
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